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The multiconfiguration Dirac-Fock method allows to calculate the state of relativistic electrons 
in atoms or molecules. This method has been known for a long time to provide certain wrong 
predictions in the nonrelativistic limit. We study in full mathematical details the nonlinear model 
obtained in the nonrelativistic limit for Se-like atoms. We show that the method with sp + pd 
configurations in the J = 1 sector leads to a symmetry breaking phenomenon in the sense that 
the ground state is never an eigenvector of or 5^. We thereby complement and clarify some 
previous studies. 
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Introduction 

Simulations of relativistic systems of Atomic and Molecular Physics and Chemistry are now 
of widespread use and the need for reliable methods is stronger than ever [24l [M]. The 
difficulties of doing relativistic simulations are however probably largely underestimated. 
One the the most common method is the so-called muHiconfiguration Dirac-Fock (or Dirac- 
Hartree-Fock) theory. This method has been known for a long time to provide certain wrong 
predictions in the nonrelativistic limit [311 [23l [131 [HI [HI [H] ■ The purpose of the present 
article is to explain in full mathematical details the encountered difficulties. 

In atomic relativistic calculations, one usually imposes the total angular momentum 

= (L+S)'^ whereas in nonrelativistic calculations both and S"^ are imposed. It has been 
observed by Kim et al in '15 that a certain multiconfiguration ground state of the symmetry 
space J = 1 for i?e-likc (Beryllium-like) atoms, was converging in the nonrelativistic limit to 
a state which was not an eigenfunction of and L^. This led to erroneous values of certain 
transition probabilities like spin-forbidden ones. 

In the multiconfiguration methods, the wavefunction is taken to be a linear combination 
of certain configurations. Both the linear coefficients and the orbitals in the configurations 
are variational parameters, leading to a highly nonlinear problem (with respect to variations 
of the orbitals). Each configuration is itself a sum of Slater determinants whose coefficients 
are fixed such that the configuration belongs to a chosen symmetry subspace ( J = 1 in the 
case of the example studied in [15]). 

It was noticed in [14] that the "error" is slowly disappearing when the number of deter- 
minants is increased. This suggests that the phenomenon is purely nonlinear, and that it 
has nothing to do with the nonrelativistic limit procedure in itself: it is the model obtained 
in the limit which does not fulfill the usual symmetry properties of nonrelativistic Quan- 
tum Chemistry or Physics models. In the limit, one obtains a nonlinear model for which 
all the configurations have J — 1 but they do not necessarily have a fixed orbital angular 
momentum or a fixed total spin S^. In those cases, the calculated ground state is not an 
eigenfunction of or S^. 

In a linear model, any nondegenerate ground state automatically has the symmetry of 
the system but in a nonlinear model there could be a symmetry breaking phenomenon: 
although the system has a certain symmetry, the ground state does not necessarily possess 
this symmetry. One then obtains several minimizers and it is only the set of all ground states 
which is invariant under the symmetry group. 

In this paper we study in detail the nonlinear model obtained in the nonrelativistic limit 
for i?e-like atoms in the symmetry J — 1, when only s, p and d shells are considerej^, 
following [IS]. Such an atom has four electrons which can only be in the following subshells: 
lsi/2, 2si/2j 2pi/2! 2p3/2, 3(^3/2 and 3c?5/2- Each subshell has a certain degeneracy but it 



That is, we only consider the lowest "ungerade" states in the symmetry J = 1. 
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is only described by one radial orbital. The distinction between the 2pi/2 and the 2^3/2 
subshells is an artefact of the nonrelativistic limit. In nonrelativistic theories, the same 
radial orbital is used for 2pi/2 and 2p^/2- The same holds for 8^3/2 and 3^5/2- The fact that 
the radial orbitals of 2pi^2 and 2^3/2 (or 8^3/2 and 8^5/2) are allowed to be different is, 
loosely speaking, similar to enlarging the variational set, which leads to symmetry breaking 
as will be explained below. 

In MCDF theory one considers all the possible configurations of the symmetry J = 1 
which can be constructed upon these states for four electrons, and one optimizes both the 
radial orbitals of the subshells and the variational coefficients in front of the configurations. 
Among all these possible wavefunctions, only few of them are eigenfunctions of both and 
S'^. In the present case, there are three possible symmetries: ^Pi, '^Pi and ^Di (correspond- 
ing to specific eigenvalues of and S"^, see Section [2] for a precise definition). In Theorem 
12.11 we give necessary and sufficient conditions on the radial orbitals and the coefficients for 
a wavefunction to be in one of these three symmetries: the radial orbitals of 2pi/2 and 2p3/2 
(and of 8^3/2 and 3^5/2) must be the same and certain relations must hold between the 
configuration coefficients. 

The issue discussed in [15l[TTl[T4] is whether 'the' ground stateH of the symmetry J = 1 
is an eigenfunction of and S*^. In other words, does it belong to one of the previous 
symmetries? To address this question, we look at the ground states obtained by imposing 
each of the symmetries ^Pi, '^Pi and "^Di, and ask ourselves whether these states can be 
stationary points and/or local minima of the full model where only ,/ 1 is imposed, or 
not. Our results fTheorems 12.21 [^75] and are summarized in the following table: 





is a stationary point, but is not a local minimum. 


'Pi 


is a stationary point, but is not a local minimum. 


'Pi 


has lower energy than ^Di and 'Pi, but is not a stationary point. 



This shows in particular that enlarging the variational set by only fixing J = 1 instead of 
both and leads to symmetry breaking for the wavefunction: the obtained new ground 
state is never an eigenvector of or S^. It is however essential that all s, p and d orbitals 
are included in the model (similar results are expected when more shells are added, by the 
same arguments as the ones presented in the paper). As we will see in Theorem 11.11 there 
is no symmetry breaking for the ground state if only s and p orbitals are considered. 

As noticed first in [15], the question whether the constrained ground states of the sym- 
metries ^Pi, ^Pi and ^Di are stationary points of the full model or not, is related to certain 
properties of the occupation numbers of the orbitals. This is explained in detail in Remarks 
O and ED below. 

The appearance of symmetry breaking is a well-known phenomenon in atomic multicon- 
figuration methods; it was encountered and rigorously examined in [2]. Similar issues occur 
in molecular calculations with regard to spacial symmetry, see, e.g., [5j [71 [22l |6] . In this pa- 
per we do not propose any practical solution to this phenomenon. Lowdin who emphasized 
this issue in a famous discussion [H] (after stability results 251 SI HSl 133 11] in Hartree-Fock 
theory) called it a "symmetry dilemma" . Our impression is that it is inherent to the way 
calculations are currently done. For the model studied in this paper, the only reasonable 
solution is probably to increase the number of determinants in order to decrease the effect 



^We put 'the' in quotation marks to emphasize that there is no uniqueness. 
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of nonlinearities. 

We have tried to make our work accessible to both the Mathematics, Quantum Physics 
and Chemistry communities. In particular, we will state and prove some well-known results 
(like for instance a special case of the first Hund's rule [16]) for the convenience of the math- 
ematical reader. We hope that our work will help in improving the exchanges between the 
different communities. On the one hand mathematicians should find the material allowing 
her/him to understand the models and the typical problems encountered in nonlinear quan- 
tum computations. On the other hand we hope physicists and chemists will appreciate our 
rigorous clarification of the phenomenon they have to deal with. 

Notation 

We recall that the angular momentum reads L = a; x (— iV), that = (Tfc/2 where Ufc are 
the well known Pauli matrices, and that J = L-|-S. For an A^-body system, we still denote by 
L = X {~i^Xk ) the (vector- valued) angular momentum operator of the N particles. 

A state will be denoted as "^^^^Lj when it is an eigenfunction of with eigenvalue S{S+1), 
of with eigenvalue L{L + 1) (with the identification P, D, F, ... for L = 1, 2, 3, ...), and 
of with eigenvalue J( J -I- 1). We will use the notation when it is an eigenfunction 

of and (with the same eigenvalues as before), but not necessarily an eigenfunction of 
J2. For more details, we refer to [TTl [29l [27] . 

1. Model vi^ith sp configurations only 

We consider Be- like atoms, i.e. atoms with 4 electrons, in the symmetry J = 1. We start 
with a simplified multiconfiguration method employing only s and p shells. This means the 
4-body wavefunction takes the form [TS] 

^' = a$'(ls2/2 2si/2 2pi/2)(i?0, i?2) + fo*'(lSi/2 2Si/2 2p3/2)(i?0, i?3) 

where a and b are configuration- mixing coefficients and i?o, i?2 and are the radial 
functions of, respectively, the shells lsi/2, '2si/2, 2px/2 a-nd 2p3/2- The two configurations 
$'(ls^^2 2s 1/2 2pi/2) and ^'{lsf^2 2s 1/2 2^3/2) are some linear combination of Slater deter- 
minants made upon the corresponding 4-component shells, with fixed coefficients chosen 
such that J = 1. We will not write the detailed form of the configurations here, but later 
we will give the precise expression of their nonrelativistic limits (see (II. 5p and (|1.6p ). The 
Hamiltonian to be used is the Dirac 4-body Coulomb operator which reads 

where — ca. ■ p + mc^f5 is the Dirac operator [5S]. 

Most atomic MCDF packages aim at calculating a certain critical point of the energy 
(i?o, ^1, a, h) I—)- (vP, H"^), under the following constraints: 

/•oo /"OO /"OO /"OO 

/ \RQ{r)\^ r^dr = / \Ri{r)\^ dr ^ / |i?2(r)|^ r^dr = / \Rz{r)\^ dr ^ I , (1.1) 
Jo Jo Jo Jo 



poo 

/ Ro(r)Ri{r) r^dr = 0, 



(1.2) 



Failure of Multiconfiguration Methods 



5 



a^ + b^ = l. (1.3) 

As the energy is not bounded below due to the negative spectrum of the Dirac operator, 
it is a priori not at ah obvious which critical point has to be considered and calculated 
numerically. Its existence in the infinite dimensional setting is also not clear at all. However, 
using the methods of O E], one can prove that well-chosen critical points of this energy 
converge as c — >■ oo to critical points of a certain nonrelativistic model which we will now 
describe in detail. 

The variational set of 4-body wavefunctions obtained in the nonrelativistic limit contains 
all the functions having the form 

^ = a$(ls2/2 2si/2 2pi/2)(i?0, Rl,R2) + fo$(ls?/2 2si/2 2p3/2)(i?0, Ri,R3) (1.4) 

but this time the two configurations are 2-component functions, i.e. they only depend on 
the spin variable. Their relation with the usual nonrelativistic configurations are given by 
the following formula (see, e.g., |3] page 294) 

1 /2 
$(ls2/2 2si/2 2pi/2)(i?o,i?i,i?2) := -^^Psp{Ro, Ri, R2) + ^^Psp{Ro, Ri, R2), (1.5) 



/2 1 
$(ls2/2 2si/2 2p3/2)(i?o, Ri, R3) ^ ^PspiRo,RuR3) + ^ 'P«p(-Ro, -Ra). (1-6) 

Here ^Psp and ^Pgp are some nonrelativistic configurations chosen such that 

(52 _ _^ 2fc+lp^^ ^ 0^ (^2 _ 2) 2k+lp^^ ^ ^j2 _ 2) 2k+lp^^ ^ 

for k = 0, 1, and which are made only of s and p orbitals. We should probably rather use 
the notation '^'^^^Pi,sp but we refrain to do so, for the sake of simplicity. The form for the 
functions ^Psp and '^Pgp (with = 1) is provided in Appendix A for the convenience of the 
reader. For many of our arguments, we will not need this explicit form. 
The associated minimization principle reads 

EspiJ=l):= inf (^-,7?*) (1.7) 

Ro;Ri.R2;R3,a,b 
satisfying ]l.ll >. il.2l and ]1.3| | 

where ^ takes the form (|1.4p and H is the nonrelativistic Hamiltonian 

• 2 \x^\l ^ \x:,-x,\' 

Note that the Hamiltonian H is real, hence each eigenfunction can be chosen to be real. 
For this reason, when passing to the nonlinear case, although in principle H should act on 
complex functions, we will restrict ourselves to real mixing coefficients a, h and real- valued 
radial orbitals R^. The same is done in most Quantum Chemistry or Physics packages. The 
extension of our arguments to complex functions does not present any difficulty. 

The Hamiltonian iJ commutes with both the total orbital angular momentum L and the 
total spin S. For this reason, in a usual nonrelativistic multiconfiguration method, one always 
restricts the search to ground states of a certain symmetry class. The only configurations 
made of s and p orbitals satisfying J = 1 are the ones corresponding to L = 1 and S — 
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(singlet) or = 1 (triplet), and which appear in ()1.5p and (|1.6|) . Therefore, we will compare 
(|1.7|) with the usual nonrelativistic methods described by the following variational problems: 

EspCPi):^ inf {'PspiRo, Ri, H 'P.piRo, Ri, R2)) , (1.8) 

Ro ,Ri , R2 

Sri r ir^fyi Ti_y, 1 11. Il l and l |1.2| | 

Esp^Pi):^ inf {^PspiRo,Ri,R2).H^P,piRo,Ri,R2)). (1.9) 

satisfying ]1.1D and ]1.2t 

Note that there is no configuration-mixing coefficient to optimize in the above minimization 
problems. For this reason, (|1.8p and (ll.9[) should indeed be called Hartree-Fock methods. 

Our wavefunction (|1.4p is always an eigenfunction corresponding to L = 1 but it is not 
necessarily an eigenfunction of S^. However, taking R2 = R3 and choosing a = (resp. 
a = —b^/2), we see that our variational set of functions of the form (|1.4p indeed contains all 
possible functions '^Psp{Ro, Ri, R2) (resp. ^Psp{Ro, Ri, R2))- Hence we deduce that 



Eap{J = 1) < min 

{EspCPi),EspCPi)}- 



(1.10) 



The specific case of (the first) Hund's rule proven below in Theorem 11.21 tells us that indeed 
Espi'Pi) < Esp{^Pi), see Corollarv ll.il In principle, however, the inequality in (|1.10p could 
be strict in which case the minimizer would not be an eigenfunction of S^, but instead a 
linear combination of S" = and S = I states. 

We will see that this problem indeed does not occur for sp mixing, as expressed by the 

Theorem 1.1 (The nonrelativistic limit for sp is correct). We have 

Esp{J = 1) = EspCPi) < Esp{ iPi). 

Additionaly , any ground state for Esp{J = 1) satisfies a = eb^/2, i?3 = eR2 for some 
e = ±1. 

The rest of this section is devoted to the (simple) proof of the above theorem. As we will 
need it in the following, we start by proving a special case of the well-known (first) Hund's 
rule (for an excellent discussion of Hund's rules, we refer to |16j . where a result similar to 
the following one is proved). 

Theorem 1.2 (Hund's rule for singlet/triplet states). For K > 2, M > I, let 

U {fi}^i be an orthonormal system o/L^(M^,C) anci 



M / K \ M / K \ 

j=i \i=i / j=i \i=i J 

where f^x,a) = /(x)5^(cr) and (cj) G C*^ \ {0}. Let 

2(K+1) 

H := (^«)/2)»+ V{xi-Xj) 

i=l l<i<j<2(K+l) 

be a Hamiltonian where h is a self-adjoint operator on L^(R'^,C) and V is a positive real 
function. Then, i/^'i and ^2 belong to the form domain of H , one has 

(*l+*2,i/(*l+*2)) < (*1 -*2,-ff(*l -*2)). (1.11) 



^Throughout the paper, we use the convention (/i A ■ ■ ■ A /]v)(^i, ■■■,^n) = (Nl) det(fi(xj)). 
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Remark 1.1. Let us emphasize that we do not impose any spacial symmetry on the func- 
tions fi and gj. 

Proof. We have 

(^-i - 1'2,-ff(*l - *2)> - (*1 + ^'2,i?(*l + ^'2)) = -45R(*i,iJ*2) = -43fi(*i,V*2> 

where V is the interaction (two-body) potential involving the function V. In the last equality 
we have used that each Slater determinant appearing in always contains two functions 
orthogonal with all the functions in any of the Slater determinants of ^2, which implies that 
the one-body term vanishes. Calculating the two-body term one gets 



(*i,V*2> = -// V{x-y) 



M 

^c, g2j^i{x)g2j{y) 



2 



dx dy 



and the result follows. □ 

Corollary 1.1 (Hund's rule for sp mixing). We have Esp{^Pi) < Esp{^Pi). 

Proof. Using the methods of proof of [HI [20l [TOl [18], one can see that there exists 
{Rq,Ri,R2) minimizing Esp{^Pi). We now choose and fix these functions. Using formula 
(jA.5| given in Appendix A and Property (jA.8| . we see that 

( ^Psp{Ro, Ri, R2) , H{ '^Psp{Ro, Ri,R2))) 

= i(S'-(3p2.sp(i?0,i?l,i?2)) , ff5-(3p2..p(i?0,i?l,i?2))) (1.12) 

where the state ^P2,sp is defined in (jA.l[) . Indeed we have precisely 

S-{^P2,sp{Ro,RuR2)) ^s^Ro)AsHRo)A(^s^Ri)Api{R2) + sHRi)A pI {R2 ) 
and, see (|A.4|) . 

^Psp{Ro,Ri,R2) = -^s^(i?o) A s^(i?o) A (^s^(i?i) Apt(^2)-s^(i?i)Apl(i?2)^, 

the notation being that of Appendix A. The result follows from (|1.12p and Theorem 11.21 □ 

We now give the proof of Theorem 11.11 

Proof. Using formulas (|1.5p . (|1.6I) . (IA.4P and (jA.5[) . we deduce that any trial wavefunction 
for (|1.7p can be written: 

* = ^ ^Psp{Ra,Ri,-aR2 + bV2R3) + ^Psp{Ra, aV2i?2 + bRa). 
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Since H commutes with S, the scalar product between the above two eigenfunctions of 
(corresponding to different eigenvalues) vanishes, and we get 

H^) = l( ^PspiRo, Ri,-aR2 + bV2R3),H ^Psp{Ro, -ai?2 + bV2R3) 



I ( ^Psp (i?o , Ri , aV2R2 + bR3),H ^Psp (Rq , Ri , aV2R2 + bi?3 ) 



> 



> II aR2 + bV2R3\\ i^^^(ip^) _ Esp^Pi)} + EspCPi) > E^p^Pi) (1.13) 



where 1 1 • 1 1 denotes the (r^ dr) norm and we have used that 



iV2i?2 + bR. 



3(a2+62) =3. 



By Hund's rule fCoroUarv ll.ip . we have Esp{ ^Pi) — Esp{ '^Pi) > 0, hence we get that Esp{J = 
1) > Esp{ "^Pi). Therefore there must be equality in p.l3l) and it holds ai?2 = bV^R^. Taking 
the square of the previous relation and using that /(i?2)^ — J {Rs)^, we prove the result. □ 



2. Model with sp and pd configurations 

In the previous section we have seen that the nonrelativistic limit of our model with only sp 
configurations was "correct" . We now study in detail the model which was considered in |15| . 
The idea is to add configurations by considering d shells. The nonrelativistic wavefunction 
now takes the form: 

* = a$(ls2/2 2si/2 2pi/2)(i?o, Ri,R2) + 6*(ls?/2 2si/2 2p3/2)(i?0, Ri,R3) 

+ C$(ls2/2 2pi/2 3d3/2)(i?0, ^2,^4) + rf*(ls?/2 2^3/2 Sdy2)iRo, ^3, ^^4) 

+ e$(ls?/2 2p3/2 3rf5/2)(Po, P3, ^5). (2.1) 

Here functions Ri are accounting for the radial part of each shell orbital, which are normal- 
ized like in ()1.2p . Only Rq and i?i have to be orthogonal. As before one can first consider 
the relativistic model with 4-component wavefunctions and pass to the nonrelativistic limit 
c — > 00. One obtains the above form (|2.ip of the wavefunction. The first two functions of 
(j2.ip have already been defined in (jl.Sp and (|1.6p . The other three functions are given by 
(see again [3] page 294) 

$(ls2/2 2pi/2 3d3/2){Ro,R2,Ri) -.^ ^%diRo,R2,R4) - ^%diR0,R2,R4) 

+ ^^Dpd{Ro,R2,R4), (2.2) 



$(ls^/2 2p3/2 3d3/2)(i?0,i?3,i?4) 



^PpdiRoi Rs, Ri) H — ^ "^PpdiRoi R3, Ri) 



15 

-—^^Dpd{Ro,R3,Ri), (2.3) 
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$(ls^/2 2p3/2 3d5/2)(^o, ^3, ^5) ^Ppd{Ro, R3,Rz) + ^Ppd{Ra, R3, R5) 

V5 vlO 



^ ^Dpd{Ro,R3,R5)- (2.4) 



The above nonrelativistic configurations satisfy, for fc = 0, 1 

(52 - A;(A: + 1)) '''+%d = {L' 2) ''^+%d = { j' 2) 2'=+iPp, = 0, 

{S^ - 2) ^Dpd = {L^ - 6) ""Dpd = ( J2 - 2) ^Dpd = 0. 

Formulas for these nonrelativistic functions of the pd shells with Jz — I are given in Appendix 
A. 

2.1. Eigenf unctions of and 

Among functions of the form ()2.ip , we will be interested in the ones which are eigenfunctions 
of and L^, i.e. the ones which have the symmetry which is imposed in nonrelativistic 
calculations. Note that, contrarily to the sp mixing studied in the previous section, our 
wavefunction is a priori not even an eigenfunction of . 

We will write \1/ G ^Pi when ^' is a linear combination of configurations ^Pgp and ^Ppd- 
We use similar notations for ^ G ^Pi and ^ £ ^Di. The following result will be crucial in 
our analysis: 

Theorem 2.1 (Eigenvectors of and 5^ of the form (|2.ip ). Let a normalized 
wavefunction of the form (j2.ip . 

(1) We have ^ G '^Pi if and only if there exists e, e' = ±1 such that R2 — e R3, R4 — e' R5, 
a ~ e\f2h, 3d = ie'e and 4c = —e\/5d. In this case 

* = ^(^\/2 + -^^ ^PspiRo,Ri,R2)~cV6%diRo,R2,R4), (2.5) 

(2) We have ^I' £ ^Pi if and only if there exists e,e' — ±1 such that R2 = sR^, R4 — e' R^ 
a\f2 — —e b, 3d = —e'e and c = —ey/bd. In this case 

^^f = ~aV3^PspiRQ, Pi, P2) + c V3 %d{RQ, P2, P4). (2.6) 

( 3) We have ^ G '^Di if and only if there exist e, e' , e" = ±1 such that a = b = 0, R2 — ePa, 
P4 = e' R5, c = e"/V^, d = J^Jhee" , e -^1/lOee'e". In this case 

^ = s"^Dpd{Ro,R2,R4). (2.7) 

Proof. Using formulas (I2.ip - (|2.4[) . we get (for the sake of clarity, we omit to mention Rq 
which appears in all configurations) 

* - 4^ ^Psp{Ri,aV2R2 + bRs) + ( %d{-cV5R2 + idR^, P4) + %d{3eR3, P5) 
V o V 30 ^ 

+ ^Psp{Ri,-aR2 + &V2P3) + ^ ( %d{cVT0R2 - dV2R3, P4) + %d{3eV2R3, R5) 
v3 v30 ^ 

' ^DpdicVT5R2 + 2%/3dP3, P4) + ^P'pd(~eV3P3, R5)) ■ (2.8) 
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Hence, using the orthogonality properties of the different configurations, we see that 'J G "^Pi 
if and only if 

-aR2 + bV2 i?3 =0, 
cV5R2 + (3e - d) i?3 = 3e(l - (P4, R5)) R3 , 
cV5R2 + (2d - e) i?3 = -e(l - (i?4, i?5)) , 

e(i?5-(i?4,i?5)i?4) =0. 

The last equation tells us that either e = or i?4 e'R^ with e' = ±1. If e = then the 
second and third equations imply that c — d — and only the first equation remains. If 
R4 — R5, then the system reduces to 

-a i?2 + bV2 Rs^O, 
cV5R2 + (3e - d) i?3 = , 
c\/5i?2 + (2d-e)i?3 =0, 
R4 — R5 . 

The second and third equations then imply that 3d — Ae. The rest follows from the normal- 
ization of i?2 and R3. The proof is similar for e' = — 1, and for ^Pi and '^Di states. □ 

In view of the above result, we now introduce the nonrelativistic ground state energies 
with sp + pd mixing 

E,p+,,CP,) -.^ inf (2.9) 

* oi the form Il2.5t 

E,p+pdCPi) -.^ inf (2.10) 

* of the form ^Toi 

Epai^Di):^ inf (2.11) 

* of the form ITTl 

Our goal is to compare these nonrelativistic energies with the one obtained in the nonrela- 
tivistic limit: 

Esp+pdiJ ^ I) := inf (2.12) 

* of the form | |2T1 

By definition, we of course have 



Esp+pd{J — 1) < min | £^sp+pd ( '^^'i ) , Esp+pd{ ^Pi) , Esp+pd{'Di)^- 



(2.13) 



The phenomenon which was observed by Kim et al. in [15) was precisely that 'the' ground 
state for Esp+pd{J — 1) (i-e. the nonrelativistic limit of 'the' MCDF ground state) was not 
an eigenfunction of L and S, hence it was not a solution of any of the problems Esp+pd{ 'Pi), 
Esp+pdi, ^Pi) or Esp+pd{ '^Di). This means that there must be a strict inequality < in (j2.13p . 
This relaxation phenomenon is itself the reason for the deficiency of the nonrelativistic limit 
of MCDF theory. It is a typical nonlinear phenomenon. 

We will now study with more details the three minimizers for ()2.10p . (j2.1ip and (|2.9p . 
We will in particular give some conditions under which ()2.13|) is a strict inequality, and we 
will check these conditions numerically in Section |3l 
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2.2. Study of the ^Pi state 

In this section, we prove that 'the' ^Pi state minimizing Esp+pd{^Pi) is never a ground 
state for J = 1, although it is always a stationary point of the associated energy functional. 

Theorem 2.2 (The ^Pi state). Let ^' G ^Pi be a function of the form jM]) minimizing 
Esp+pdCPi) defined in (j2.10|) . 

(^ij T/ie associated mixing coefficients and orhitals (a, ...,e, Rq, R<^) satisfying the relations 
of Theorem \2.1\ (Assertion 2) provide a stationnary point of the total energy functional 

(a,...,e,i?o,-,i?5) ^ (*,if*> 



where ^ takes the form (|2.1|) . 
(2) However, ^ is never a local minimum of this functional. 

Proof. We first show that is a critical point of the total energy functional. For simplicity 
we assume that c ^ 0, the proof being the same otherwise. We also assume for simplicity 
that e = e' = 1, hence R2 = i?3 and i?4 = R^. We have to consider both variations with 
respect to mixing coefficients, and to orbitals. The vanishing of the variation with respect to 
the mixing coefficients (a, e) is a simple consequence of the fact that there is no overlap 
between states belonging to different symmetry spaces, 

for any radial functions in the corresponding states. 

We now turn to the variation with respect to orbitals. By cxtremality of the singlet 
function among ^Pi states and due to the constraints on Pq: we hav^ 



dR2 dR3 
and 



e span(Po, Pi) for i = 0, 1 , (2.14) 
R3=R2 G span(P2) (2.15) 



ap 



\R.^=B.i e span(P4). (2.16) 



dRi dR5 

Using e — 1, a \/2 — —6, 3d — —e and c = —V^d as given by Theorem 12.11 we see that 
at \1/ we have 

^pr"2^'^-'=«- (^-^^^ 

and 

^pr"3^'^==«*' ^^-^^^ 

From (PT^ - (PT71) and (pi^ - ipiB . wc find 

af(*) 



dR, 



e span(Pi) for i = 2,..., 5 (2.19) 



*The derivatives appearing below are the coordinates of the gradient of the energy with respect to the scalar 
product of I/^([0, 00), dr). 
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which ends the proof of the criticahty of 4'. 

Remark 2.1. The exceptional relation (|2.17p holds true because, as was noticed first in 
\1 5| /. for the singlet state the ratio between the occupation numbers of the Pi/2 function i?2 
and the occupation number of the p^/2 function in the ^Pgp state, is the same as the 
corresponding ratio for the ^Ppd state. Let us explain this with more details. Considering a 
variation SR G (i?2 = ^3)^; we find that the variations of the total energy functional are, 
using p.Sp . 

= 2//?* , 'P„(-Ri,M) + '-^ 'PM^R. Ri)\\ 

f|> (M, . 2(h* , + ^^^^^ «.))). 

When the matrix 

( -a cVW \ 

V3 V30 1 
b\/2 -d\/2+3eV2 I 

is not invertible, its columns are colinear and we have ^§^^{5R) — k ^^^^(SR) for some k 
(k — 1/2 in our case). Hence 

^(SR) + ^(6R) - ^ ~dR^^^'^^ = °' 
dR2 OR, ^'''^-'^\'-m^SR)^0, 



dR. 



as we want. This holds true when 

V3 _ Vso 



bV2 -dV2+3eV2 ' 
V3 V30 

Using the relations between a,...,e provided by Theorem \2.1l we see that the above equality 
reduces to (when a,c ^ 0) 

1 10 
3. — 30. 

2 20 

3 30 

which is precisely the ratio between the occupation numbers as mentioned before. The argu- 
ment is the same for Rj^ — R^. 

We now turn to the proof that ^' is never a local minimum, which will simply follow 
from Hund's rule. Note that we can write the relations between a, 6, e as 

^) = -aVs {^^ , j - cV3 I 

where Usp and Upd are the Condon-Shortley unitary matrices j3j 
1 / /2 1 \ 1 ^ 



VVI5 -V5 VTOy 
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Now we define the following new mixing coefficients 



"'^ = -aVsu* il] , U' 1 = cVsu: 



and note that by construction (a', e') is orthogonal to (a, e) for the scalar product of 
R^. Also we have 

^' ^{a',b',c',d',e',Ro,Ri,R2 = i?3,^4 = i?5) 

i.e. it takes exactly the same form as ^E* but with triplet states instead of singlet states. Now 
we vary the mixing coefficients as follows Vl — e) + t(a' , e'), which results into a 

variation for the wavefunction of the form Vl — t-^'if + W' . Calculating the energy of this 
new wavefunction we find 



Note that there is no first order term since is a stationary state as shown before (or simply 
because 4* and 5*' belong to different symmetry spaces, hence {H'i', ^P') ~ 0). 

Now we claim that f (^'') < £{'^), which will clearly imply that cannot be a local 
minimum. By (IA.8I) and (jA.lsp in Appendix A, we know that £'(5'') = where 

= -^S- ''P2,sp{R0,Rl,R2) + ¥2,prf(i?0,i?2,-R4) 

which is the simple triplet state taking the same form as but with the adequate signs 
reversed. The last step is to apply Theorem 11.21 with the following functions: /i = s{Ro), 
91 =Pi(i?2), 92 = (cV20 + a2)-i/2(c/(2V5)do(i?4) +as(i?i)), 33 =P-i(i?2), .94 = d2{Ri), 
95 = Po{R2) and ga — c?i(i?4). This ends the proof of Theorem 12.21 □ 

Remark 2.2. When R2 — ±i?3 and R4 — R5, the problem consisting of varying only the 
mixing coefficient essentially reduces to that of finding the eigenvalues of the nonrelativistic 
Hamiltonian matrix, i. e. the matrix of H in the space spanned by the 5 configurations built 
upon the orbitals: 

^Psp{Rq, Rl, R2), ^Ppd{Roy 1^2, R4), ^^'spl-Roi -R2), ^Ppd{Ro, R2, R4), ^Dpci{Ro, R2, Ra)- 
This 5x5 matrix is block diagonal and its eigenvalues are: 

XiCPi)^ inf £{a''P,p{Ro,Ri,R2) + l3%d{Ro,R2,R4)) 

A2('=Pi)= sup £{a''Psp{Ro,Ri,R2) + P%d{Ro,R2,R4)) 

^i^Di) = (^H^Dpd{Ro,R2,R4), ^Dpd{RQ,R2,R4))- 

What we have used in this second part is that Ai(^Pi) < Ai(^Pi), by Hund's rule (indeed, 
we even have that the 2x2 Hamiltonian matrix of ^Pi states is smaller than the one of ^Pi 
states, in the sense of quadratic forms). However although it is expected that in many cases 
Ai(^Pi) < X{^Di), there is no general rule: this may depend on the orbitals Ro^.-.^Rq. 
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2.3. Study of the ^Di state 

In this section, we prove that 'the' ^Di state minimizing £'pd( ^Di) is also always a stationary 
point of the associated energy functional and we give a condition implying that it is not a 
local minimum. 

Theorem 2.3 (The ^Di state). Let 5* G ^Di be a function of the form (|2.7p minimizing 
Epdi^Di) defined in ([2lT|) . 

(1) The associated mixing coefficients and orbitals (a, ...,e,Ro, ...,i?5) given by Theorem \2.1\ 
provide a stationary point of the total energy functional 

(a,...,e,i?o,-,i?5) ^ (2.20) 

where takes the form (|2.1|) . 

(2) If moreover 

inf £ {a^PspiRo, Ri, R2) + P^PpdiRo, R2, R4)) < Epd^Di), 

fl.l6{flo}^, l|flllL2(,2^,, = l 

(2.21) 

where Rq, R2, R4 are the radial functions of , then the stationary state 5* is not a 
local minimum of the functional (j2.20p . 

Remark 2.3. Our assumption (|2.2ip exactly means that \{^Di) = Epdi^Di) is not the 
lowest eigenvalue of the Hamiltonian matrix as explained above in Remark \2.2\ . It is not clear 
how to prove (|2.2ip rigorously. In Section\3[ we will verify it numerically on the approximated 
solutions provided by the ATSP HF and MCHF program by Froese-Fischer Jl^ . 

Proof. Wc do not give all the details of the proof which is very similar to that of the ^Pi 
state. The fact that is a stationary point is seen exactly as in Theorem 12.21 

The proof that 5" is not a local minimum is also very similar to that of the ^Pi state, 
with the difference that we do not have a general Hund's rule for ^Di states, hence we need 
to require condition ()2.21|) . which is expected to be true in many cases. This time we have 
(assuming again e — e' — e" — 1 for simplicity) 

The result is then obtained by arguing as before with, this time, 

where a and (3 are chosen to minimize the left side of (|2.2ip . □ 

2.4. Study of the ^Pi state 

In this section, we give a simple condition implying that 'the' '^Pi state minimizing 
Esp+pdi^Pi) is not a stationary point of the total energy functional. 
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Theorem 2.4 (The '^Pi state). Let 'J G "^Pi 6e a function of the form (|2.5P minimiz- 
ing Egp+pdi^Pi) defined in (j2.9p . and denote by {a, ...,e, Rq, R5) the associated mixing 
coefficients and orbitals satisfying the relations of Theorem \2.1\ 
If there exists 6R E (Ps)^ such that 

c{H^, 3Pj,d(Po, <5P, P4)) ^0 or a{H^, ^PspiRo, Ri,SR)) ^ 0, (2.22) 

then 4' does not provide a stationary point of the total energy functional 

(a,...,e,Po,...,P5) ^ 

Remark 2.4. Condition (|2.22p is wer?/ intuitive. It indeed implies that 

Esp+p'd{^Pi) < Esp+pd{^Pi) , (2.23) 

where 

Esp+p'dCPi)= inf £(a¥,p(Po,Pi,P2)+/3^Ppd(Po,P3,i?4)). 

(i?0,---,-R4) 

satisfying il.ll and ( II. 2t 

T/ie relation between (|2.22p and (j2.23p was already noticed in 

Note that when a,c ^ 0, one can prove that (ifvP, '^Ppd{Ro,SR,R4)^ ^ or that 
(iJ^P, ■^P<,p(Po, Pi, (5P)) 7^ for some SR and these two conditions are indeed equivalent. In 
Section[^ the condition (P^, ^Psp(Po, Pi, (^P)) ^ is verified numerically. 

We now give the 

Proof. As before, it can easily be seen that our state ^> is indeed a stationary state with 
respect to variations of the mixing coefficients only. The non-criticality will come from the 
variations of the orbitals, as suggested by (|2.22p . As before we assume for simplicity that 
e = e' = 1, hence P2 = P3 and P4 = P5 for the ground state ^. 

By extremality of the triplet function ^ among '^Pi states, we have similarly as before 

a£(^') 



aP2 dR3 
and 

a£(*) a£(*) 



espan(Po,Pi) for i = 0, 1 , (2.24) 
R3=R2 e span(P2) (2.25) 



OR, 



fls=7?.4 G span(P4). (2.26) 



9P4 aPg 

This time we find using the relations of Theorem 12.11 that 

dRi ^ ^pT 

hence only variations with respect to P2 and P3 remain to be considered. The main point 
is that there is a priori no relation between ^q^^ and ■ More precisely, let us consider 
a variation ^P e (P2 = Pa)^. We have 

^(^P) = 21h^ , (^-^ ^Psp{R,,5R) + ^ 'PpdidR, P4)^ ^ 
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and 

= 2(h^, (^^3p,^(R,,SR) + ^^%,{SR,R^)^'^, 

where we have used the relations of Theorem 12. II in the last line. The main difference with 
the ^Pi state is now that the matrix 



(2.27) 



is invertible. Hence we get that ^' is stationary with respect to variations of i?2 and R3 
independently if and only if 

ySR e iR2 = Rs)^ , a{H^, ''PspiRi,6R))=c{H^, ^Ppd{SR, R^)) = 0. 

This clearly leads to a contradiction when (I2.22p holds true. □ 

Remark 2.5. The fact that the matrix ()2.27p is invertible can be interpreted in saying that 
the ratio between the occupation number of the pi/2 function R2 and the one of the P3/2 
function R3 in the sp configuration is not the same as the one of the pd configuration U5f . 



2.5. Conclusion 

In the previous sections we have studied the states of the nonrelativistic symmetries '^Pi, 
^Pi and ■^£'1. As a consequence of our results we obtain the 

Corollary 2.1 (Occurrence of symmetry breaking for sp + pd). Assume that (12.211) 
holds for a ^Di ground state and that condition p.22p holds for a "^Pi ground state. Then 
we have 

Esp+pd{J = 1) < niin |i?sp+pd( ^Pi) , Esp+pd{ ^Pi) , Esp+pd(, '^^i)|'- (2.28) 

Additionally a ground state ^ for Esp-\-pd{J = 1) is never an eigenfunction of and neither 
ofS\ 

Proof. The strict inequality p.28p is an obvious consequence of the previous results. It 
implies that a ground state for the minimization problem Esp+pd{J = 1) cannot be a com- 
mon eigenfunction of and S^. What remains to be proven is that it cannot even be an 
eigenfunction of or of separately. This means that any ground state must have a 
nonvanishing projection in each of the symmetries ^Pi, '^Pi and ^Di. 

Assume for instance that 5* is an eigenfunction of L^. It cannot be a '^Di state by ()2.28p . 
hence one must have ^ E ^Pi + "^Pi. Using ()2.8p . one obtains the conditions 

f cVl5R2 + (2dV3 - eV3(P4, P5)) ^3 ^ 0, 
{ e(P5-(i?4,P5)i?4) -0. 
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If e = we get R2 — eRz and c\/b — ~2ed for some e e {±1}- Hence 



^ ^ ^^^^ ^Ps,{R,,R2) + ^%4R2,R,) 



— ^ ^ ^Psp{Ri, R2) H — ^Ppd{R2, R4) 



By Hund's rule we get 



+ eb ^ 



V3 



'Psp{Ri,R2) 



V3 



Ppd{R2, R4) 



Esp+pd ( ^Pl ) 



+ 



^Psp{Ri, R2) H — ^ ^Ppd{R2, R4) 



Esp+pdi^Pi) > Esp+pdi^Pi) 



which contradicts ()2.28p . The argument is the same if R4 ~ ±i?5. Hence we have shown that 
cannot be an eigenfunction of L^. The proof that ^ cannot be an eigenfunction of S*^ is 
very similar. □ 



3. Numerical verification of (P^^ and ^^I^ 

In this section we verify numerically the two assumptions made in the previous section. We 
use the package ATSP of Froese-Fischer [T2] to get approximations of the nonrelativistic 
states in the different symmetry spaces. The program uses 220 discretization points on a 
logarithmic grid and a finite difference method. We have then verified conditions (|2.2ip and 
((2^ with the help of Mathematica. 



3.1. Verification of ((2:22)) 

First we start with the MCHF calculation for the configuration ^Psp+pd for the Beryllium 
atom. This yields constants a,c and functions Rq, Ri, R2, R4, which are numerical approx- 
imations of the real ones. We can also run the program for the simpler Hartree-Fock cases 
'^Psp, ^Ppd and we get the following total energies (in Hartree). 

S,p(3Pi) ~ -14.5115, 
E.p+pd^'Pi) ~ -14.5166. 

Also we obtain for the ^Pgp+pd configuration the following numerical values for a and c: 

a - 0.9951963, c ~ -0.0978997 . 

The plots of the radial parts of the orbitals Rq, Ri, R2 — R3 and R4 — R5 (corresponding 
respectively to the shells Is, 2s, 2p and 3d) are displayed in Figure [1] below. 

Now, we are going to verify (|2.22p by showing the existence of a function SR, orthogonal 
to i?2, such that 



(i?*, ''Psp{Ro,Ri,5R))^0. 



(3.1) 
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Fig. 1. Plots of Ro (top left), Ri (top right), R2 (bottom left) and R4 (bottom right) for the optimal 
^Pi{sp + pd) configuration of the Beryllium atom, obtained with the package ATSP. 



As explained in detail in Appendix B, we can calculate the exact expression of (j3.ll) in terms 
of the radial functions only. We obtain 



00 poo 2 j-2 

/ (2|i?oGs)p + R2{t)6R{t) ^ . dsdt 

Jo max{s,t} 

°° r°° s^t^ minis t} 

/ Ro{s)5R{s)R,{t)R,{t)- ^pJ-dsdt 

Jo "J max \s, r| 

a / Ri{s)6R{s)R2 (t) R, (t) ds dt 

Jo Jo "J max \s, t\ 

cV2 / Ri{s)R2is)SRit)R^{t) ^ .^tL^dsdt. 

Jo Jo 3max^{s,i} 



Finding a function 6R, orthogonal to R2, for which the above formula is away from 
is less convincing than arguing as follows. The above formula is linear with respect to SR 
hence it can be written 



/•oo 

{H^, ^Psp{Ro,Ri,SR)) = / F{r)SR{r)dr 

Jo 
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with 

/ \ R (r) 

F{r) ■.= a{~-R'^{r)~rR',{r)+R2ir)-4rR2{r))+a / {2\Rois)\' + \Ri{s)\') \ . ds 
\ ^ /Jo max-[s,rj- 

poo 2 2 • r 1 poo 2 2 • r 1 

Then we note that the vanishmg of F{r)SR{r) dr for all SR orthonormal to R2, is indeed 
equivalent to the existence of a Lagrange multiplier A such that F = Ai?2. Hence it suffices 
to show that F{r)/R2{r) is not constant in order to verify (I2.22p . 
We compute numerically F/R2 and get the following graph: 




Fig. 2. Plot of F{r)/R2{r) for r between 1 and 5. 



This clearly shows that the function F/R2 is not constant. Hence (|2.22l) holds true. 
3.2. Verification of (|23T|) 

Let us now pass to the numerical verification of p.2ip . In order to do so, we again run 
the Froese-Fischer ATSP program jT2! for the HF calculation of the "^Dpd configuration for 
the Beryllium atom. This provides us with radial functions i?o, R2 and i?4 for the orbitals 
Is, 2p, 3d. For the 2s orbital, we construct a function i?i by taking the one obtained in the 
calculation for the ^Psp+pd configuration and projecting it on the orthogonal to the space 
generated by Rq . We get the following numbers for the total energies (in Hartree): 

Epdi^'Di) ~ -14.1889, 

£ {^PspiRo, Ri, R2)) ^ , 

thus yielding a numerical verification of condition (|2.2ip . 

Appendix A. Nonrelativistic configurations in the sector J — 1 and — 1 

For the convenience of the reader, we quickly explain how to construct the nonrelativistic 
configurations used in the text. In the whole appendix we use the notation 

s^R){x,a) = Ri\x\)6ria), pl,{R){x,a) = R{\x\)Yr f ^r(^), 
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-1 





1 


-1 


3P2 = 4 si Apii 


2^2 si A pi 


2 si A pi 





iPl = 2^sT Apij- si ApL^j 
2^sT Apij + si Apli^ 


lPl = y2(^sTApi-siApJ^ 
y2(^sT A pi + si ApJ^ 


iPl = si ^ 

si Apj; 


pj — si A pj 
+ si Apj 


1 


2 si Apli 


sT A pj < 


\- 1 

3P2 = 

L- 


s- 

: si Apj 



Jz = 1 Jz = 2 



Table 1. Derivation of all the common eigenvectors of , 5^, Lz and Sz, starting from the ^P2,sp state 
having all the highest possible quantum numbers, at the bottom right of the table. For the sake of clarity, 
we have not indicated the closed shell Sq{Ro) A Sg(-Ro) A • ■ ■ which appears in front of all the configurations. 
We have also refrained from mentioning the radial functions R\ and R2. 



dl,{R){x,o) = R{\x\)Yi^^[^^^5r{a) 

for T G , i} and m ~ — ^, Here F/" are the usual eigenfunctions of and Lz (orbital 
angular momentum), normalized in the Hilbert space i2(S'2), such that 

and which can be expressed in terms of Legendre polynomials [3] . 
sp configurations 

We start with sp configurations. Table [1] summarizes the different possible states which are 
eigenfunctions of both and S"^, classified with respect to their and Sz- The usual 
method is to start from the state having the highest possible quantum numbers (S" = 1, 
= 1, L = 1, = 1, J = 2 and Jz — 2) which is located at the lower right corner: 

^P2,sp{Ro,Ri,R2) -.^ s^Ro) A s^Ro) A (^s^Ri)Apl{R2)y (A.l) 

Next we apply the lowering operators — Lx — iLy and = Sx — iSy to get two 
states having the same 5 = 1 and L = 1 but without a precise J: 

5-(^P2,sp(i?o,i?i,i?2)) =s^(i?o) As^(i?o)A (^.s^(i?i) Apt(i?2) + s^(i?i)Apl(i?2)), (A.2) 
L-{^P2.sp{Ro,Ri,R2)) :=y2/(i?o)A5^(i?o)A (^/(i?i)ApJ(i?2)). (A.3) 
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A usual antisymmetrization of (jA.2|) gives our (normalized) singlet state having S = 0, 
hence J = L = 1: 

^Psp{Ro,Ri,R2) = -^s^Ro) A sHRq) A (^s^Ri)Api{R2)-sHRi)Apl{R2)y (A.4) 

In the subspace — 1, we can construct a '^P2 state by applying the lowering operator 

J- = L- + S-: 

J {^P2,sp{Ro, Rl, R2)) — L~ (^^P2.sp{Roi Rl, R2)) + (^^P2,sp{Ro, Ri, R2)) ■ 

We know that the subspace = 1 is of dimension 3. Hence we deduce by orthogonality 
that our (normalized) '^Pi state with J = 1 and = 1 is 

'^Psp{Roj Ri, R2) = ^ ^-^P2,sp{Ro, Rit R2) ■ (A-5) 

Note that 

L+L-(3p2,.p(i?o,i?i,i?2)) =2(3p2,sp(i?o,i?i,i?2)), (A.6) 

5+5-(¥2,,p(i?o,i?l,i?2)) =2(¥2,,p(i?o,i?l,i?2)), (A.7) 

which explains why 1/2 is the right normalization in (jA.sp . From (|A.6p and (|A.7p we see 
that for any observable A commuting with L and 5* (for instance A = H, our nonrelativistic 
Hamiltonian) , then we have for all a, b and all Rq, Ri, R2 and R'q, R'l, R'2 

(aL- +65-)(3p2,.p(i?o,i?i,i?2)) , A (aL- + bS-){''P2,spiR'o, R'l, R'2)) 

{C? + 1?)^- ( 3P2.sp(i?0, i?2)) , A L-CP2.sp{K,R[,R2)) 

= (a2+62)/5-(3p2,,p(Po,-Ri,i?2)), A5-(3p2,,p(P^,p;,P^))\ 

\ / L2(R3)4 

= 2{a^ + h')hp2,sp{R^^Ri,R2). A^P2.sp{R'o,R'i,R'2)) ■ (A.8) 

\ / L2(R3)4 

pd configurations 

We now switch to the calculation of the pd configurations. 

Table [2] shows the dimension of the different common eigenspaces of and Sz for pd 
configurations. As before we may start from the ^Fj^ state having the highest quantum 
numbers and derive our states in the J = 1, = 1 sector by applying successively the 
lowering operators P", 5*^ and . 

The final result is as follows. The ^Pi state in the = 1 sector reads 

%d{Ro,Ri,R2) = ^s^(Po) As^(Po) A (^V6pl,{Ri) A d^^{R2) - VSpliRi) A d\{R2) 

+ pI{Ri) a 4(P2) - \/6pii(Pi) A 4(P2) + %/3pi(Pi) A 4(P2) -p\{Ri) A ^(Pa) 

(A.9) 
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_3 


_2 


_1 


n 

yj 
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2 




-1 


1 


2 


3 
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3 


2 / 

3 r/ 3^ 


1 / 

/ 3p 





2 


4 


6 


6 


/ 6 y 


'Z 4 / 

3^ 3 n/ 


/ 2 / 

li?/f 3p 


















1 


1 


2 


3 


/ 3 / 


/ i / 







































J, = l = 2 = 3 J., = 4 

Table 2. Dimensions of the different common eigenspaces of Lz and Sz for configurations. In eacfi cell, 
we have indicated the type of state spanning the corresponding space. 



The ^Pi state reads 



— 

Ppd{Ro, Ri, R2) ^ 7; P2,pd{Ro, Ri, R2) 



(A.IO) 



where 



'P2,pd(i?o, i?2) = s^Ro) A s^Ro) A (^pl.iRi) A 4(i?2) - -^pliRi) A 4(i?2 

VvlO vlO 



V3 



+ —p[{R,)Adl,{R2)] (A.ll) 



is the unique configuration having = 2. Note that, like in the sp case, one has 

L+L-{^P2,pd{RQ,Rl,R2)) = 2('V2,prf(i?o,i?l,i?2)), (A.12) 



S+S-{^P2,pd{Ro,Ri,R2)) ^2{^P2,pdiRo,Rl,R2))■ 



iA.13) 



Therefore we deduce that (jA.SP stih holds for linear combinations of the triplets sp and pd 



Failure of Multiconfiguration Methods 23 

and that for any observable A commuting with L and S, 

a^Psp{Ro, Ri, R2) + b^Ppd{Ro, R3, R4) , A (a! ^Psp{Ro, R[,R^) + b' ^Ppd{Ro, -R3, -R4)) 
1 



= ^{^aS-{''P2,sp{Ro,Rl,R2))+bS-{''P2,pd{Ro,R3,R4)) 

A {a'S-{^P2,sp{R'o,R[,R2))+b'S-{^P2,pd{Ro,R3,R4))) 

(A.14) 

= 3P2,.p(i?o, i?2) + b^P2,pd{Ro, Rs, R4) , 

A {a' ^P2,sp{K, R[, R2) + b' ^P2,pd{R'o, R3, K)) ) ■ (A.15) 

/ L2(R3)4 

Lastly the '^Di state is obtained through the more complicated formula 

^Dpd{Ro, Ri,R2) = ^ {L-L- - 3L-S- + 65-5") ^D3,pd{Ro, -Ri, -R2) (A.16) 
where ^Ds^pd{Ro, Ri, R2) is the unique "^D^ state having 3, given by 
''Ds,pd{Ro,Ri,R2) = sHRo)AsHRo)A(^^pl{Ri)AdUR2)-^pl{Ri)Adl{R2)y (A.17) 

Appendix B. Energy expressions as functions of the radial components 

In this appendix, wc provide the formulas of the energy written in terms of the radial 
components of the orbitals. We will first need the 

Lemma B.l. Assume that we have 6 mutually orthogonal functions, /i, . . . , /e- Then, 

^/iA/2A/3A/4, (^73^) -^i^-^^A/sA/e^ = <^/3A/4, ^^/sA/e^, 

(B.18) 



/1A/2A/3A/4, I ]_ I /i A /2 A /3 A /6 



/1A/4, ^^/lA/6\ + //2A/4, ^^/2A/6\ + //3A/4, ^-^^3 A , 

\x-y\ / \ \x-y\ / \ \x-y\ / 

(B.19) 

and 



/1A/2A/3A/4, Y I 1 I /l A /2 A /3 A /4 



i<i<j<4 ^ \ y\ / 
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The proof of the above lemma is based on long but straightforward computations. We 
use the above result to compute some quantities that are needed in order to perform the 
numerical computations of Section [3] 

We start with the computation of {^Psp{Ro, Ri, R2) , H ^Psp{Ro, Ri, R2)) which is 
needed to verify (lOTj) . By (|X8)) we have 

( ^Psp{Ro, Ri, R2), H ^Psp{Ro, i?2)) = ( ^P2.,sp{Ro, Ri,R2) , H ^P2,sp{Ro, Ri,R2)) ■ 
A simple calculation shows that 

( '^P2,sp{Ro,Rl,R2) , Pt ^P2,sp{Ro, Rl,R2)) 



^ {2\so{Ro)ix)\^ + \soiRi){x)\^ + \piiR2){x)\^) dx (B.22) 

+ hp2,sp{R0,Rl,R2), i 1.1 I ''P2,sp{Ro,Rl,R2))- (B.23) 

The well-known properties of the harmonic spherical functions show that 



(IB^TI) + (Ib:22]) = ^ k'i?o(r) 



2 2 



R2{r)\' 



2\ 



dr 



4r {2RQ{rf + R^{rf + R2{rY) dr. 
On the other hand, using Lemma FB .11 the last integral ()B.23|) is equal to 
mm ^11 ko(i?o)(:c)n5o(i?o)(y)P + (2|so(i?o)(x)p + bi(i?2)(x)P) |so(i?i)(y)|' dxdy 



2\s^{R^){x)\'\pi{R2){y)\' dxdy 



{so{Ro){x)so{Ro)iy)+pi{R2){x)pi{R2)iy)) so{Ri){x)so{Ri){y) dx dy 



So (i?o) (a;) So {Rq ) {y)pi {R2){x)pi (i?2 ) (y) dx dy . 

Using the well known formulae that can be found for instance in Slater's book [27^ (Section 
13-3 and Appendix 20a), this can be rewritten as 



(Ro{s)^Ro{t)^ + (2i?o(s)' + R2{sf) Ri{tf + 2Ro{sfR 
Ro{s)Ro{t)Ri{s)Ri{t) 



{tf 

§2 ds dt 



s 



2 ds dt 



max{s, i} 



^0 



(i?2(s)(i?2W {Ri{s)Ri{t) + Ro{s)Ro{t))) 



max {s. i\ 

s^t^ vtmY{s,t \ dsdt 



3 max^ {s, t\ 
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We then go on to calculate the expression of (i/^*, '^Psp{Ro,Ri,SR)), appearmg in 
(|2.22p . As a corollary of Lemma IbH of (|A.5|) . and of (|A.8p . we obtain 

(^PspiRo,Ri,R2), f 5Z I ^ I I ^Psp{Ro,Ri,6R)) 

= \P2,sp{Ro, Ri, R2) , ( ^ I 1 I ^P2,sp{Ro,Rl,SR)\ 

\ Vi<»<j<4 i^^^^^ V / 

{2\so{Ro){x)\^ + \so{Ri){x)\^) pi{R2){y)pi{SR){y)dxdy 
So {Ro){x) Pi (SR) (x) pi (i?2 ) ) So {Ro){y) dx dy 
so{Ri){x) pi{dR){x) pi{R2){y) so(i?i)(?/) dxdy . 



R3> 

Using again the formulae in Slater's book [37] (Section 13-3 and Appendix 20a), we finally 
deduce that 

hPsp{R^,Ri,R2)A I I ^Psp{Ro,Ri,SR)\ 

r<x> poo 2 j-2 

= / / {2\R„{s)\^ + \Ri{s)\^) R2{t)6R{t) ' dsdt 
Jo Jo max {s,E| 

/ Ro{s)6R{s)R2{t)Ro{t) ^^^dsdt 

Jo Jo "J max \s, t| 

- r r Ri{s)SR{s)R,{t)R,{t) "^^^''^^ dsdt. 
Jo Jo 3max^{s,i} 

On the other hand, by (jA.15|) . and using the same methods as above, we find 

{ %d{Ro, i?2, i?4) , — r '^'*p(^o, Ri,sR) ) 
\ Vi<»<j<4i^' "^^7 / 

= ( ^P2,pd{Ro, R2, R4) , I 1 T 1 ^P2,sp{Ro, Ri,SR)i 



, X 7 X -7 I 

^l<i<J<4 ' ■" 

V2 / / Ri{s) R2{s) 5R{t) R^{t) \ "? .1 dt . 

o max \S ^ tj 



Jo 
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